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I. INTRODUCTION 


An outer-loop control system synthesis procedure is currently under- 
going development. A description of the approach is given in Ref. 1. 

The outer-loop synthesis procedure is being applied to a 3-D R/NAV control 
design problem for the NASA ATOPS B-737 research aircraft. 

The outer-loop design approach requires an operational inner-loop 
control system before the outer- loop feedback gains can be determined. 

The inner-loop control system employed is a derivative the Boeing control 
system documented in Refs. 2-7. The advantage of using the Boeing design 
is that is has been extensively analyzed and flight tested. 

To facilitate the design process, a linear model of the inner- loop 
control system is required. This technical note presents the linear 
and nonlinear models of the Boeing inner-loop control system. Validation 
of the linear model is determined by comparing linear and nonlinear 
simulations. 

There are two approaches which can be used in representing the 
inner-loop control system. The easiest approach, and the method used in 
this report, is to model the inner- loop control system in continuous time 
and combine the resulting equations with the linear models of the aircraft 
and actuators, to obtain 



The states are the n-vector, the outer-loop controls are the m- 

vector j , and the disturbances are the d-vector , «, , 

“*cl cl cl 

and are constant matrices of appropriate dimension. The discrete 
model for the outer— loop design becomes 

^K,k+l “ ^CL^^L.k *^CL^^L,k 


( 2 ) 


where 


^CL = ® 




^CL = e ds e ds 


(3) 

(4) 


The sampling time, At» can be chosen as desired. 

The second and more accurate approach Is to model the plant and 
actuator dynamics in continuous-time as before but model the inner 
control system in discrete time. The assumptions used in digitally 
implementing the inner- loop control system (trapezoidal integration, 
measurement and computation delays, Tustin's transformation, sampling 
time = 0.05 sec) are accounted for in the discrete model. A discrete 
representation of the plant is combined with the discrete inner-loop 
control law system to produce Eq. 2. If the outer-loop control system 
is to operate at a slower rate (At « 0.1 sec), the control design becomes 
multirate. Techniques for handling multirate multivariable control 
designs have recently been determined in Ref. 8 but are outside the scope 
of this work. The second approach represents future schemes which may be 
useful if more accurate evaluation techniques become necessary. 
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The report is organized into six chapters. Chapter two presents 
nonlinear block diagrams of the inner-loop control laws and the corre- 
sponding linear time-invariant equations of the inner-loop models. 
Chapter three details the steps used in linearizing the inner-loop 
control laws. Chapter four combines the linear inner-loop models into 
the continuous-time plant representation shown in Eq. 1. The properties 
of proportional plus integral feedback of acceleration error in the 
inner-loop control laws are also discussed in Chapter four. Chapter 
five presents simulations comparing linear model time histories with 
nonlinear model time histories. The report is summarized in Chapter 
six. Appendix A compares nonlinear ACSL simulations with a recently 
developed ATOPS B-737 FORTRAN simulation. 
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II. INNER-LOOP MODELS 


A. AUTOTHROTTLE 

The autothrottle inner-loop is the most con^licated of the control 
designs. Two error signals are computed and compared to threshold 
values. Control feedback paths are switched on and off as the error 
signals cross the threshold values. A wlndshear estimator Is used to 
produce feedback for wlndshear compensation. The versine of the roll 
angle Is fed back to advance the throttle just as the airplane rolls 
into a turn. The roll angle signal is washed out to prevent throttle 
advance in a steady turn. A block diagram of the autothrottle inner- 
loop control system is shown in Fig. 1. The gains in the block diagram 
are defined in Table 1. 

A detailed explanation of the wlndshear estimator is given in Ref. 

3. Basically, true airspeed, TAS and longitudinal acceleration, 

are employed in a complementary filter to produce V where 

Cr 


V = -V 

CF W 3^ 

^10 


(5) 


s + 1 


V,, is the wind velocity making V a washed-out estimate of wind. V „ 
W Cr Cl? 

is next passed through a wind turbulence filter to filter out high 


frequency components in 

• ^19 1 ^ 

V = Bound (V - V „) 

CF s -k ^ CF CF^ 


( 6 ) 


The negative wind shear estimate,VcF, is subtracted from Vgs to form 
the acceleration feedback signal. The acceleration feedback signal is 
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subtracted from the acceleration command, x , to form the command error. 

c 

The integral of the command error becomes one of the feedback paths to 
the Incremental throttle command 

If the engine pressure ratio (EPR) approaches the maximum safe 
value f or EPR (MXEPR), special logic, (switch B as shown in Fig. 1) switches 
the EPR error signal to input the throttle integrator. The effect is to 
cause a decrease in the incremental throttle command. Switch A causes 
a throttle down error command to immediately decrease the rate of change 
of the throttle Integrator. 

A linear model for autothrottle is as follows 


4 

Av 
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(7) 


0.0 
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( 7 ) 


A6t(, =[kg 0.0 1.0 O.o] 


AV 


CF 


+ [-kg 0.0 0.0] 


AV, 


AC, 


^'^TAS'^'^’^TAS 

A<H-Av. 


L^'^WO. 


+ [0.0] [Ax^] 

The state A^^ is the Integral of perturbation command error. The state 
A<()^q is the perturbation washed out versine roll correction factor. The 
acceleration measurement, AV^g, is the perturbation rate of change of 
ground speed. AV^^g is the measured perturbation true airspeed. A^ is 
the measured perturbation roll angle and Ax^ is the perturbation outer- 
loop command. In straight and level flight, the versine roll correction 
state can be eliminated. The variable, Y, in Eq. 7 is an option to 
account for switch A as follows 


e^^„ = MXEPR - 0.1 - EPR 

jCiirK O 


z = 


0.0 


®EPR ^ ^3 

< ^EPR < ^3 


^EPR < 


OPTION 1, A > 0, Y = Z 

OPTION 2, A < 0, Y = 1.0 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


For a given trim flight conditions, outer-loop gains could be determined 

first for OPTION 1 then for OPTION 2. If noticable changes are apparent, 

the effect of A may have to be accommodated in the outer-loop control 

law. The nonlinear bounds caused by the gains k, , k, , , k and k„ are 

D li o 
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neglected v/hen the linear model is constructed in Eq. 7, An approach 
using describing functions shown later in Eq, 18 could be used to model 
nonlinear bounds in the linear plant representation. The variables 

and Av. are measurement noise. 

TAS cp 

B. VERTICAL PATH 

A block diagram of the vertical path inner-loop for the elevator 
is shown in Fig. 2. The gains are given in Table 2. Neither the 
elevator inner-loop, nor the autothrottle inner-loop employ pitch 
attitude feedback for stability during path tracking. Reference 3 
explains that replacing a pitch attitude command system with a 
vertical acceleration command system allows the aircraft to weather 
cock vertically upon encountering vertical gusts and shears so that 
glideslope beam tracking performance can be enhanced. Other types 
of control design procedures can produce an elevator control system 
that weather cocks, employs pitch attitude feedback, and has good 
glideslope tracking performance as discussed in Ref. 9. 

The elevator inner-loop filters the vertical acceleration measure- 
ment, forms the acceleration command error, then feeds the error signal 
in proportional-integral form to the elevator. The integrator output 
is position limited. Washed-out pitch rate is fed back for improved 
stability. The vertical acceleration command is position limited to 
insure passenger comfort. The feedback signal to elevator is multiplied 
by a gain that decreases with increasing calibrated airspeed. The 
linear model for elevator is 
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( 13 ) 


" *^CAS^*'27 *^25**^22 ^25^ 


"^%J0 

/\ 

Ah 




Aq + Av 
Ah + Av;; 


*^CAS^^25*’*^22^ 


(14) 


The perturbation states are the integral of command error, A^- , filtered 

h 

/N 

vertical acceleration, Ah, and washed out pitch rate, A^^. The measure- 
ments are pitch rate, A^, and INS (Invertial Naviation System) vertical 
acceleration. Ah, The outer-loop control input is the vertical acceler- 
ation command, Ah^, The limit on Ah^ could be handled similiar to Eq. 18, 

C. HORIZONTAL PATH 


The horizontal path inner-loop feeds back roll rate for stability 
augmentation and the error between the outer-loop roll angle command and 
roll angle. A block diagram of the control law is shown in Fig. 3. The 
gain values are defined in Table 3. The feedback signal which forms the 
aileron actuator command is multiplied by the gain k^ which is a function 
of calibrated airspeed, CAS, as shown in Fig. 3. The roll command signal 
is limited. 
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A linear model of the horizontal path inner~loop which accounts 
for the roll command rate limit is 


A^‘ 


-1.0/k^2 


1.0/k32 0.0 


Ap + Av 

P 

/s 

S= 


-h 




Me 


0.0 -k3gj A(|>^^ 


0.0 0.0 


A4> + Av . 

9J 








Ap 1 + [0.0 j"Ap + Av - 

A$cJ Lm + 


+ [0.0] [A(})^] 


( 15 ) 


(16) 


The perturbation states are filtered roll rate, Ap, and filtered roll 
command, A(|)p. The measurementas are roll rate, Ap, and roll angle A({)* 

V* 

The outer-loop control input is the roll command, A({)^. 

Using describing functions, if A(j)^ changes abruptly then 

A = A(|)^ - A$^ = 4>c - 


1 1.0 

1.0 

IT 


The advantage of allowing J to be variable is that outer-loop guidance 
gains can be designed for different values of A. 


a| < k 


37 




4 2 


+ 2 sxn 


'|A|) 


(18) 


A > k 


37 


D. RUDDER 

The Boeing inner-loop control system for rudder is the yaw damper 
shown in Fig. 4. The gains are shown in Table 4. Body axis yaw rate 
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is filtered to suppress measurement noise. The filtered yaw rate is 
multiplied by a gain which is a function of CAS as shown in Table 5- 
to decrease the gain as airspeed increases. The filtered yaw 
rate signal is washed out for turn coordination then position limited 
to reduce control authority. Outer-loop commands directly actuate the 
rudder surface. The Boeing outer— loop control design commands rudder 
for decrab during landings. 

The linear model of the inner-loop is 
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Ar, 


WO 


Ar 


+ [0.0] [Ar + Av^] + [1.0] [Au^] 


(19) 

( 20 ) 


The position limits in Fig. 4 could similarily be incorporated in the 
linear model using Eq. 18 if this is deemed desirable. The perturbation 
states in the linear model are filtered yaw rate, Ar, and washed out yaw 
rate, Ar^^. The measurement is yaw rate, Ar. The outer-loop control 
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III. LINEAR MODEL ANALYSIS 


A. LINEARIZATION OF FEEDBACK ELEMENTS 

The Boeing inner-loop system uses the derivative of ground speed, 

• •• 

true airspeeds and vertical acceleration, h as feedback 

Go lAb 

elements. Linear analyisis requires that these elements be linearized 
and expressed in terms of the perturbed states and controls of the air- 
craft. This section derives the perturbation relationships. 

Figure 5 shows the relationships between the accelerations of the 
vehicle in a local-level north pointing frame and the along track and 
cross track accelerations. The north pointing frame is denoted as the 
geographic coordinate system. In the figure, V is the ground speed, a 
is the vehicle acceleration in the horizontal plane, and 5 is the ground 
track angle. 

From the figure it follows that 


= 2 . 


tan - ^ , 

X 


„ ..2 . . 2 .% 

Vg = (x + y ) 


(21a, b) 


The along track, a.„^, and cross track, a_„_, accelerations are related 
AliN. GIK 

to X and y through the transfomation. 


ATK 


CTK 


[■ cos? sln?1 rxl 
tsin? cos?J LyJ 


( 22 ) 


Expressing cos? and sin? in terms of x and y the following occurs 


ATK 


• •• , 9 .• 

_ XX + yy 


(23) 
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CTK 


- xy - yx 


(24) 


Note that taking the derivative of V in Eq. 21b, shows that V„ and 

G G ATK 

in Eq. 23 are equal as required. Perturbing Eq. 24 produces 


AVgH 


X - V_ X v_ y - y 
Go o Go o Go Go 


Go 


Go 


^ "] 


Ax 

Ay 

Az 


Ax 

Ay 

Az 


(25) 


which is rewritten in matrix form as 


= H Ax + H.. Ax 

G X — X — 


(26) 


Let H ((j),6,if)) be the transformation from body axes to geographic 

G 

axes where (j), 6, and ijj are the Euler angles (platform axes and geographic 
axes are assumed to coincide) . It follows that 
T 




h; 


B 


• „G 


(27) 

(28) 


where 


Vg = [u V w] 


(29) 


and 


„G • „G --G 

^ = «B ^B »B “b ^B 


(30) 
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The vector uyl represents the body axis angular rates in radians, (a flat 
B 

nonrotating earth is assumed) 


\:^<n 

li 

[p q 

r] 


o5_ is the matrix 
j5 

representation of the 

by 

0 

,1 


~G 

r 

0 

1 



-q 

p 0 


Perturbing Eq. 28 produces 

o 

- <*o> 

4 h 

O Q 

where 

’l.O 

0.0 

-slnQ 

o 

S = 

o 

0.0 

cosd) 

slnd) COS0 
o o 


0.0 

-sin4 

o 

cosd) COS0 
0 


(31) 


(32) 


(33) 


(34) 


and 


AVp = [A(j) A0 Ai|)] 


(35) 


Perturbing Eq. 30 produces 


AE = AVg - Kg AVg + AVg 

o o o o 

+ I-B AVb - Hb ^'^SSb 

o o o 


(36) 


A expression for Av„ can be determined from the linear aerodynamic model 

.jg 
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of the aircraft. 


Ax = A Ax + B Au + E Aw 
- - - -g 


( 37 ) 


Extracting the equation for Av„ from Eq. 37 produces 




[AAA 3 
vv vw vv 


Av, 


Aw. 




+ B Au + E Aw 
vu — vw 


(38) 


The states. Aw are gust disturbances. Substituting Eq. 37 into Eq. 36 
results in 


AH * (Hg w®) AVg + (H® A^ - Vg ) Aw® 

o o o o o 


«V*o> “b H S ®B 

o o o o o 


+ HI B Au + E Aw 
B vu — B vw — c 
o o 


(39) 


Substituting Eq. 33 and Eq. 39 into Eq. 36 determines the desired equation 
for Vg in terms of the aircraft perturbation states and controls. 


AV_ = H Av„ + H w_ + H Av„ + D Au + D Aw 
G w — B vw— B w ~B vu — vw — 


(40) 


where 


H = H.H„ + H..H, 
w X B X 1 
o 

H = H® A + H® w^ 
1 B vv B B 
o o 


H = H..H- 
vw X 2 


(41) 

(42) 

(43) 
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( 44 ) 


o o 


"w - %“3 

o o 

H- = [ (x -X ) - X ] H^ L„ + H? F 
3 o o B B ^ W 

o o 


D = H..H, 
vu X 4 


H, = H“ B 
4 B vu 
o 


D = H..H. 
vw X 5 


H. = E 
5 B vw 
o 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


A similiar expression for Ah is obtained by noting that 


Ah = H^Ax 


(51) 


where 


= [0 0 - 1 ] 

Substituting Eq. 39 into Eq. 51 determines the desired perturbation 
expression for Ah in terms of the perturbation aircraft states and 
controls. 


(52) 


Ah = + \H 3 AVg + H^H^Au + Hj^H 5 Aw (53) 

True airspeed is the velocity of the vehicle relative to the 
atmosphere, 

''tab ■ <“a ''a *A>** <54) 
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u* 


V. 


w. 


V. = V - w- - w 
—A — B G 


( 55 ) 


(56) 


The vector ^ represents the 3-axls steady state wind above the earth's 

— iii 

surface and is not modeled in this analysis. Frcan Eq. 54, the perturbed 
value for is 


AV, 


TAS 


- “VA 


(57) 


where 


»VA 


J ^Ao 

V 

TAS 

u. r 


V 


Ao 

/ 

TAS 


^Ao ~ | 

\as J 


(58) 


Substituting the perturbed value for AV obtained from Eq. 56 produces 

A 

the final expression for 

''\as - “va - V As! 

«• 

The final expressions relating the perturbed values for AV^, Ah, 
and AV^^g and the linear aircraft model are determined by distributing 
the elements in Eqs. 59, 53, and 40 into Eq. 77 for the chosen order 
of Ajc and A^ in Eqs. 78 and 79. 


B. GUST MODEL 

The gust terms in the model are of a random nature and can be 
Modeled using the well-known Dyrden spectrum. Ref. 10. The modeling 
effort consists of using spectral factorization methods to obtain a 
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dynamical system which generates a random process having the specified 
power spectral density when driven by a white noise process. Ref, 11, 
Rotational gusts around the aircraft are ignored. The transfer functions 
for the gusts are as follows. 


u = 
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rr 

(1 + L s) 
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TAS 


V 
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(1 + L s)^ 

V 

V 

TAS 
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w 



1 + v/r^ 
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(1 + L 
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w 


TAS 


(60) 


(61) 


(62) 


The airspeed, V_.„, is defined in Eq. 54, L , L , and L are the scales 
^ ’ TAS’ ^ u v’ w 

of turbulence, and a , o , and a are the variance of the gust. The 

u V w 

scales and gust variance are shown in Table 6. 

A state-space realization of Eqs. 60 to 62 is 


Aw = 



= H Aw 
w -g 


(63) 





( 64 ) 
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1 0 0 0 0 
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in Eq. 64 is a 3-vector of independent Gaussian white noise processes 
with unit variance. 
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IV, A LINEAR MULTIVARIABLE INNER-LOOP 


MODEL USEFUL FOR OUTER-LOOP DESIGN 


The block diagrams for the inner- loops in Figs, 1, 2, 3, and 4 
are relatively straight forward and easily transformed into a nonlinear 
simulation using the ACSL programming system. Formulating a linear 
multivariable model for the inner-loop system which is useful for linear 
simulations and linear analysis is the purpose of this section. The 
acceleration feedback in the inner-loop complicates the derivation. 

Each inner- loop linear model discussed in Chapter II can be placed 
in the following form 

= ElAx^ + E2[A^j^ + IL^ 

A6g^ = E4Axg + E5[Aj^^^ + Av^^] 4* E6Au^ (69) 

The elevator inner-loop model is used as an example. The other inner- 
loop linear models representations use A, R, or T in place of E in 
Eqs. 68 and 69 for aileron, rudder, or throttle inner-loops. The states 
of the inner-loop dyanmic models shown in Chapter III are as follows 



/S 

Ah 


(70) 

a4 - 

/V 

AV^ 

''♦wo' 

(71) 

= [Ap 


(72) 


Av] 


(73) 
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The vector represents the commands to the inner-loop system 


AuJ - [Aic^ Ah^ A4>^ Au^] 


( 74 ) 


These commands become the new controls when the inner-loop control system 
closes the loop around the aircraft dynamics. The vector A^ represents 
the inner-loop measurements used for feedback. 


T 

A^ = [Ap Aq 

Ar AO AV ^3 AVj ^3 AhJ 

(75) 

T 

'"'tas 


Aa - [AVgg 

(76) 


The last three inner-loop measurements, A^, are expressed in terms of 
the aircraft states, controls, and gust in Chapter III, 


Aa = C Ax + D Au + D Aw (77) 

~a--au— aw— 

The vector Asc represents the states of aircraft, in body axes, 

= [Au Aw Aq A0 Av Ap Ar A<|> A^ Ax Ay Az] (78) 


The vector, Aji, represents the aircraft controls, 

a/ . AS^ AA^J,^ A6^I (79) 

The vector Eqs. 68 and 69 are white zero-mean Gaussian noise 

states representing the measurement noises for the sensors used in the 
inner-loops. The vector has the following covariance 


E{Av AVj ) - DIAfiONAL [O? 0?_ O. 


Ap “Aq “Ar “A* “AVgg “Av^^^ '^Ah' ’ ’iL 


The expression 0 . is the standard deviation for the measurement noise 

Ap 

of the roll rate gyro. The other variables in Eq. 78 are standard 
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deviations of the measurement noises for the pitch ratio gyro, yaw rate 
gyro, roll angle from the INS (Inertial Navigation System), inertial 
along track acceleration from the INS, airspeed sensor, and vertical 
acceleration from the INS. The actuator command vector, A^, is composed 
of elements from each of the inner-loop control systems discussed in 
Chapter II, 

“ec- “rc> 

The model for the linearized actuator dynamics is discussed in Ref. 13, 
and has the following form. 


Ah 

Ax 


F X + G A6 
u— u u — c 

^ A X + B A6 
u— u u — c 


(82) 

(83) 


In Eq. 79, A6„„ is thrust while A6_„ in Eq. 80 is the throttle command, 
ill iC 

If all the inner- loop models are combined, the result is as follows. 


+ ®IY t^^IL ^^IL^ ^ ®IC^^c 

A^ = [AXjj^ + Avjj^] + (85) 

AXil = + °IL^^ 

Table 7 shows how the matrices in Eqs. 84 to 86 are constructed. 

In the rest of the derivation, the inner-loop control system is 
closed around the aircraft perturbation dyanmics. Substituting Eq. 86 
into Eq. 85 and the result into Eq. 81, one obtains 



where 


fp T ' T * T T * 

Ax = [Ax Aw Ax Axi: ] 
— u — IL 


K 


11 


[GuHiyCil 


F G ] 
u u IL 


^21 = ^u^IC 


^31 = ^u«IY 


K# - s= G 
41 u lY IL 


Defining 


Z„ = (I - K41) 


-1 


Eq, 87 reduces to 


Au = K^Ax + K^Au + K.Av«y 
— 1 — p 2 3 — IL 


where 


^1 = "u^ir S = Khv h = "uSi 


Rewriting Eq, 86 as 


AVtt = K, Ax 4* K- Au 
-^IL 4y ^ 5y 


where 


V - [<=IL 0 “IW 01 


K. - D„ 

5y IL 


and substituting for A^ from Eq. 94 results in 


( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 
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( 99 ) 


AVtt = K,Ax + K_Au + K,Av_^ 

“^IL 4 — p 5 — c 6 — IL 

where 

S = Sy^2 

K, = K- K- 
6 5y 3 

Substituting Eq. 99 Into Eq. 83 and regrouping produces 
AXil “ ■*' ^®IY ®IY^6^^-IL ^®IC 

where 

Aj . [0 0 0 Ajj^] 

Substituting Eq. 99 into Eq. 85 produces 

A4 - 1% + ”lAlA£p + [Hiv + + [Hj^ + Hj^jlAu^ 

where 

= [0 0 0 Hjj^] 

Substituting the above into Eq. 82 one obtains 

A4 - [Aj + B^[Hj + 

where 

A2 = [0 0 0] 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 

(107) 

(108) 
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The aircraft dynamics satisfy the equation 


= AAx + BAii 

Substituting Eq. 94 into Eq. 109 produces 

Ax = [A„ + BKjAx + BK-Av„ + BK„Au 
— 3 1 ~p 3 — IL 2 — c 


(109) 


( 110 ) 


where 


Ag = [A 0 0 0] 


( 111 ) 


Combining Eqs. 110, 64, 107, and 103 determines the desired closed- loop 
model using the inner-loop control system for the linear aircraft 
dynamics 


■ M ' 


■ Aj + BKj 

•* 



Aw 


0 A 0 0 

w 


Aw 

+ 

Ax 
— u 


*2 “u"*! + 

Ax 
— u 


/iiL 


_Ai -h B^^K^ 

J 


L 



"bk^ 


! 

”o " 


+ 

0 

[ASti 

r + 

B 

w 



BuEHiy + 

L ” 

u 

0 



.»IY 



0 


BK„ 


'iC + Vs 


H 


An 


( 112 ) 


Equation 112 is the expanded version of Eq. 1 discussed in the introduction. 


A. FEATURES OF ACCELERATION FEEDBACK 


The Boeing inner-loop control system feeds back the integral of the 
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acceleration error qunatities 


G*<* 

X 



•• 

X 

c 


(113) 


anct 


= ii 


ii 

c 


(114) 


to improve stability and tracking. The properties of this type of feed- 
back can be studied by investigating a simple scalar system 

X = ax + bu (115) 

with the control law 


u = kj^(x - x^) + k^ 5 
5 = /(x - x^) 


(116) 

(117) 


The closed- loop system for the scalar plant is 


[%] 

The two rows in the closed-loop plant matrix in Eq. 118 are linearly 
dependent implying that one of the closed-loop eigenvalues is always 
zero. Two of the four zero eigenvalues in Table 9 are caused by the 
way integral feedback is used in the Boeing inner-loop control law. 

The steady-state tracking ability of the control law can be 
investigated using the model 


X 


a4-bC 


ki bk^ 

l-k^b^ l-k2b 
k, bk„ 


a+b(- 


1-k^b' l-k^b 



— 


-bki 


X 

1 

1-kj^b 


c 

-r 

-bk^ 



1-k^b 
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(119) 


X 


a 0 


X 

+ 

b 

u + 

0 



1 0 


K 


0 


-1 



. - 


- . 






when is a constant ^ an expression for u can be obtained from Bq, 116. 
The closed- loop system for Eq. 119 using u is 




ak-+k„ 




f-bk 1 

•• 


^ ii_ / f\ 


. 



X 




X 


1-bk, 


sr 

1 




1 

1 


1 0 


c 


-1 

^ ~ 




» J 


» » 


The steady state value for x is given by 


( 120 ) 


bk. 


X = 


a+bk^ ^c 


( 121 ) 


Unless a is zero, there is no value for which makes x equal to x^ in 
steady state. A step command for Ah^ in Chapter IV shows that perfect 
steady state tracking is not obtained in simulation even though integral 
feedback is employed in the Boeing inner-loop control system. 
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V. INNER-LOOP MODEL VERFICATION 


The purpose of this chapter is to compare time histories between 
the linear inner- loop closed-loop models and the nonlinear ACSL simu- 
lation model which includes nonlinear models of the inner-loops control 
systems. A number of options are available when the linear model is 
constructed. Table 8 shows the options and the recommended settings 
currently employed. 

The effect of varying the engine dynamics EPR time constant is 
investigated in the longitudinal dynamics verification. The effects 
of the lead/lag filter and spoiler aerodynamics are investigaged in the 
lateral-directional dyanmics verification. The choice of the first two 
options shown in Table 8, remain unresolved. 

A. LONGITUDINAL 

The simulation comparison has two purposes. The first purpose is 
computer coding verification. The second purpose is to identify non- 
linearities which cause a significant descrepancy between linear and 
nonlinear dynamics. 

The autothrottle has a significant nonlinearity in the upper and 

lower saturation limits for the rate of change of EPR (engine pressure 

ratio) as discussed in Ref. 12. The effect of the saturation limit is 

to decrease the EPR dynamics time constant in a manner similiar to Eq. 18. 

Unsaturated, the EPR dynamics time constant is -0.2 sec. A time constant 

of -2.0 sec is used in Ref. 12. Figure 6 shows the effect of T for 

EPR 
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-0.2, -0.5, and -2.0. The simulation comparisons in Fig. 6 are in good 
agreement. The recommended T is -0.2, i.e., no saturation effect 

UilrK 

is needed in the linear model. 

Originally there were significant mismatch between linear and non- 
linear autothrottle simulations. The mismatch was traced to a subtle 
error in the calculation of the upper saturation limit of EPR in the 
ACSL nonlinear simulation, which was subsequently corrected. The error 
in the ACSL program existed, but apparantly did not affect the results 
in Ref. 12. 

The time history comparisons for a Ah^ command is shown in Fig. 7 
and has excellent agreement. The Ah response is also shown in Fig. 7 
and as discussed in Chapter III, the steady state error is not zero 
even though the inner-loop control system feeds back the integral of 
the command error. 


The inner-loop closed-loop eigenvalues for the longitudinal system 
are shown in Table 9. The phugoid mode and the short period are stable 
and well damped. The effect of changing T__ has almost no effect on 

Ctl: K 

the eigenvalues. Residualizing the elevator actuator dynamics primarily 

affects the eigenvalue for the Ah filter by further stabilizing its 

value. The two zero eigenvalues for A^ and A^- are a feature of feeding 

X h 

back the integral of acceleration as discussed in Chapter III. The states 


associated with the gust model (Aw Aw and Aw and the wind shear 

gl g2 g3 

estimator (AV^^ and AV^) are uncontrollable and have eigenvalues which 
remain fixed for variations in the aircraft dynamics. 


B . LATERAL-DIRECTIONAL 


The largest discrepancy between linear and nonlinear models occurs 
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for the aileron inner~loop control system^ The discrepancy is caused by 
the highly nonlinear aileron/ spoiler actuator system discussed in Ref. 12. 
The linear model for the inner-loop control system and actuator is con- 
structed so that any subset of the nonlinearities can be represented in 
the model. 

Three linear models are simulated for step A(J)^ commands of 2.0 deg 
and are shown in Fig. 8. The roll command is large enough so that spoiler 
is activated when A<{)^ = 2.0 deg, but when A<|)^ returns to zero, little 
spoiler is used. The assumptions in constructing the three linear 
models are discussed in Table 10, which shows the closed-loop eigenvalues. 

The closed-loop eigenvalues indicate that the Dutch Roll mode, spiral 
mode, and roll mode are very stable. In CASE 1, the slow actuator eigen- 
value, -0.975, is deceptive since it can be shown that the eigenvalue is 
almost cancelled by a zero of nearly equal value. The aileron actuator 
in CASE 1 responds almsot immediately to the commanded aileron value as 
long as the value has a finite rate. The simulations in Ref. 12 for 
step (% infinite rate) commands in the aileron actuator can be misleading 
in this regard. 

CASE 2 models the aileron actuator as it existed before many of the 
mechanical nonlinearities were added. The fast aileron actuator mode 
forms a complex pair with the -yaw rate filter mode. The roll mode 
becomes more stable. 

CASE 3 includes the effect of spoiler. In the nonlinear model, 

spoiler is activated after the aileron actuator surface commands exceeds 

certain values (% 2.2 deg). Including the effect of spoiler when A<}) is 

c 

small produces incorrect results. Not including spoiler when A<1>^ is of 
moderate value may not produce correct results. 
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The simulations in Fig, 8 show that CASE 1 and 2 match the non- 
linear simulation for small Acf)^ commands. Using a Ac|)^ of -25 degs and 
Eq, 18, Fig. 9 shows the descrepancies that occur as A(})^ is increased. 

The effect of the saturation limit in reducing the roll response rate 
limit for large roll commands is clearly evident in Fig. 9. The roll 
command saturation limits are apparently not used in the simulations 
in Ref. 6 where -25 deg roll commands are used to test for limit cycles. 

Rudder commands of 2.0 deg are simulated in Fig. 10 for two cases. 
The rudder linear and nonlinear response match well. Continued problems 
with the linear model for the aileron/spoiler actuator system is evident 
in the aileron response in Fig. 10. Including spoiler feedback increases 
the aileron command control effectiveness. 
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VI. SUMMARY 


«9 

The linear and nonlinear inner-loop control system models to be used 
in a 3-D R/NAV outer-loop control synthesis problem are presented in this 
report. Most of the modes with fast eigenvalues are retained in the 
linear model since they do not pose a problem to the outer-loop (limited 
state feedback) control synthesis procedure. Closed-loop eigenvalues for 
the inner- loop control system discussed in the report show that all inner- 
loop complex modes are well damped and inner- loop real modes are acceptably 
stable. 

Small step commands in each of the outer-loop control variables show 
good agreement between linear and nonlinear model time histories. The 
roll inner-loop control system is Identified as the model with the most 
descrepancy primarily because of the highly nonlinear aileron/spoiler 
actuator. 
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APPENDIX A 


Recently, an alternative FORTRAN simulation of the ATOPS B-737 
aircraft was made available to ICS by NASA, As part of the verifica- 
tion of the new simulation, inner-loop step commands time histories 
computed using the ACSL simulation were similarily computed using the 
FORTRAN simulation. The descrepancies between the simulations were 
identified and partially eliminated. Figures 11 to 14 show the 
comparisons for the same step commands made in Figs. 6, 7, 8 and 10. 
The most notable difference is for the autothrottle control system in 
Fig. 11. The aileron response difference in Fig. 13 is caused by the 
fact that the FORTRAN simulation at the time the simulation was per- 
formed did not include the lead/ lag compensator and did have spoiler 
feedback. In Fig. 14, the FORTRAN simulation included the lead/lag 
compensator and spoiler feedback and the aileron responses are in 
reasonable agreement. 
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TABLE 3. AILERON/SPOILER BLOCK DIAGRAM GAIN VALUES 


GAIN 

VALUE 

'‘31 

20.0 

^32 

0.05 

^33 

1.4 

^34 

50.0 

^35 

25.0 

^36 

5.0 

^37 

4.0 

^38 

2.0 

^39 

97.66 


TABLE 4. RUDDER BLOCK DIAGRAM GAIN VALUES 



3 






TABLE 5 BREAKPOINTS FOR YAW DAMPER GAIN VERSUS AIRSPEED 


INPUT 

CAS 

(kts) 

100.00 

122.4 

150.0 

206.0 

450.0 

OUTPUT 

KYD 

1.0 

0.765 

0.61 

0.395 

0.31 


TABLE 6 SCALES AND VARIANCE FOR GUST MODELS 




ALTITUDE 
h (ft) 

0-60 

60-328 

328-1750 

>1750 

a 

u 

(ft/sec) 


16.0 

60.72h“®‘^^ 

10.0 

10.0 

o 

V 

(ft/sec) 


12.7 

26.50h"^*^® 

10.0 

10.0 

a 

W 

(ft/sec) 


9.3 

9.3 

9.3 

9.3 

L 

u 

(ft) 


560.0 

560.0 

10.85h°'^® 

1750.0 

L 

V 

(ft) 


320.0 

102. 17h^’^^ 

10.02h°'^® 

1750.0 

L 

w 

(ft) 


174.0 

12.65h°'^^ 

10.0211®'^® 

1750.0 
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OPTIONS 


RECOMMENDED VALUE 


INCLUDE SPOILER FEEDBACK 


CURRENTLY OPTIMAL 


INCLUDE LEAD/LAG DYNAMICS 


CURRENTLY OPTIMAL 


ENGINE EPR TIME CONSTANT 


-0.2 SEC 


RESIDUALIZE SPOILER ACTUATOR STATE 


INCLUDE EFFECTS OF PCU RATE LIMITS 


FALSE 


INCLUDE EFFECTS OF BACKLASH IN AILERON 
ACTUATOR MODEL 


FALSE 


RESIDUALIZE RUDDER ACTUATOR STATE 


TRUE 


RESIDUALIZE ELEVATOR ACTUATOR STATE 


FALSE 


INCLUDE EFFECT OF EPR REDUCTION IN CHOOSING 
Y IN EQ. 11 

RESIDUALIZE Ap STATE 


FALSE, Y 


FALSE 


RESIDUALIZE Ax STATE 


FALSE 


RESIDUALIZE Ar STATE 


FALSE 


INCLUDE SATURATION DESCRIBING FUNCTION 
FACTOR IN A(}) 


FALSE 


INCLUDE SATURATION DESCRIBING FUNCTION 
FACTOR IN Ax 


FALSE 


INCLUDE SATURATION DESCRIBING FUNCTION 
FACTOR IN Aii 


FALSE 
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CASE 


"^EPR 

EPR=-0.5,A6e 

ACTUATOR STATE RESIDUALIZED 

^EPR 

Ax 

0.0 

0.0 

0.0 

0.0 

Az 

0.0 

0.0 

0.0 

0.0 

AC 

X 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

PHUGOID 

MODE 

-0.'0'57±j0.03 
Wn“0 . 065 ,C=0 .878 

-0.057±j0.03 

Wn=0.065,C=0.878 

-0.057±j0.03 

Wn=n.065.C=0.878 

-0.'058±j0.03 
Wn=0 . 066 ,C=0 .878 

Aw , 



-0.138 

-0.138 

-0.138 

-0.138 


-0.149 

-0.149 

-0.149 

-0.149 


-0.149 

-0.149 

-0.149 

-0.149 

AV 

cf 

-0.20 

-0.20 

-0.20 

-0.20 

Aq 

wo 

-0.242 

-0.254 

-0.254 

-0.335 j0.155* 

'^^PR 

-5.53 

-2.11 

-2.11 

o)j^=0.369, =0.907* 

SHORT PERIOD 
MODE 

-1.77±jl.93 

0)n=2.62,C=0.675 

-1.77±jl.94 

0)n=2.63,C=0.673 

-1.78±jl.86 

Wn=2.60,C=0.693 

-1.76±ji.94 

Wn=2.62,C=0.675 

/*■ 

AVi 

-5.0 

-5.0 

-5.0 

-5.0 

<1 

-7.39 

-7.39 

-8.49 

-7.39 

A6 

e 

-23.5 

-23.5 

— 

-23.5 


* Two Reals Formed a Complex Pair 


TABLE 9 


INNER-LOOP LOGTITUDINAL CLOSED-LOOP EIGENVALUES 










CASE 

1 

2 

3 

Ay 

0.0 

0.0 

0.0 

Alp 

0.0 

0.0 

0.0 

Aw T 



-0.149 

-0.149 

-0.149 


-0.149 

-0.149 

-0.149 

Ar 

wo 

-0.273 

-0.275 

-0.397 


-0.975 

-6.09±j2.13* 

-0.975 

DUTCH 

ROLL 

-0.875±j 0.605 

-0.902±j 0.574 

-1.01±j0.24 

oj =1.06,5=0.822 
n 

0) =1.07,5=0.843 
n 

0) =1.04,5=0.973 
n 

MODE 

SPRIAL 

MODE 

-1.6 

-1.6 

-1.6 

ROLL 

MODE 

-2.73 

-3.56 

-4.58 

ROLL 

COMMAND 

-5.0 

-5.00 

-5.0 

Ar 

-5.85 

U) =6.45,5=0.943* 
n 

-10.5±j5.0* 

Ap 

-17.4 

-22.5 

0) =11.6,5=0.903* 
n 


CASE 1 

NO SPOILER AERODYNAMIC EPEECT 
IN B MATRIX, LEADLAG COMPENSATOR 
DYNAMICS INCLUDED IN A MATRIX 


CASE 2 

NO SPOILER AERODYNAMIC EFFECT 
IN B MATRIX, LEADLAG COMPENSATOR 
DYNAMICS NOT INCLUDED IN A MATRIX 


CASE 3 

SPOILER AERODYNAMIC EFFECT IN B 
MATRIX, SPOILER ACTUATOR STATE IS 
RESIDUALIZED, LEADLAG COMPENSATOR 
DYNAMICS INCLUDED IN A MATRIX 


* Reals Combined to Fom a Complex Pair 


TABLE 10 INNER-LOOP LATERAL /DIRECTIONAL CLOSED-LOOP EIGENVALUES 


•SI 




ENGINE PRESSURE RATIO 
EPR 
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FIGURE 1 AUTOTHROniE INNER-LOOP CONTROL SYSTEM BLOCK DAGRAM 
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FIGURE 2 ELEVATOR INNER-LOOP CONTROL SYSTEM DIAGRAM 
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FIGURE 4 YAW DAMPER AND RUDDER ACTUATOR MODEL D 
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FIGURE 6 (CONTINUED) SIMULATION RESPONSE FOR A Ax COMMAND OF 

2.0 FPS2. c 
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FIGURE 7 SIMULATION RESPONSE FOR A aH COMMAND OF 2.0 FPS^. 
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FIGURB.8 SIMULATION RESPONSE FOR A A<^c COMMAND OF 2.0 DEG. , 
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